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SUMMARY AND CONCLUSIONS 


Calculations have been carried out to predict unsteady-aerodynamic 
sources of wing structureborne vibration for various modes of propeller 
installation. 

Mounted either upstream of the wing or downstream of the pylon that 
supports it, a propeller at cruise must chop through a steady but nonuniform 
inflow field. For a pusher propeller the largest source of nonuniformity 
is the wake of the upstream planform; for a tractor propeller it is the much 
smoother, and therefore less severe, upstream interference of the wing at 
the propeller disk plane. Both types of inflow generate a net harmonic 
propeller thrust, collected at the hub and acting to drive the wing or 
pylon in roughly the flight direction through the engine connections. 

In the case of a tractor propeller the significant mechanism for 
generation of harmonic forces is the unsteady-aerodynamic reaction of the 
wing to the unsteady downwash induced by the convected propeller wake. 

Both configurations and mechanisms have been analyzed here for operating 
conditions and geometric parameters corresponding to the Hamilton Standard 
high-efficiency propeller. On theoretical grounds we conclude that the 
upstream wing effect will not result in a significant source of vibration. 

Our calculations in this part of the study required that we first model 
the interference wing flows of the Gu3 fstream-IT aircraft at the proposed 
propeller location. 

For a pusher propeller cutting through its pylon's wake we have predicted 
a 40-lb thrust at the fundamental (blade-passage) frequency. The model of 
wake momentum deficit for a nonlifting pylon was here based on available 
full-scale measurements for the wing of a typical airliner; it used the 
average of pairs of data points above and below the wing plane, with the 
wake .properly scaled down and applied at downstream stations where the 
propeller disk will likely be. Whether or not the 40- lb number will result 
in an unacceptable noise level depends on wing/cabin transfer functions not 

fully available yet, either theoretically or experimentally, and on whether 
the force can be isolated effectively, for example, through properly designed 
engine mounts. 
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We did not analyze fuselage interference or boundary- layer effects? 
however, for both wing- and aft-mounted propellers those flows should be 
bounded in degree of nonuniformity by the wing interference and pylon wake, 
with corresponding harmonic loads which fall between those predicted here. 

Our estimates of wing pressure distributions due to wake impingement 
also await structural/acoustic calculations for a final conclusion. 

A few general recommendations strictly based on the aerodynamics 
problems examined can be made, however: 

Non uni form Propeller Inflow as a Mechanism of Thrust Generation 

1. Having a large number of blades is a good idea, not just from the 
steady-performance viewpoint, which seeks lower individual blade tip speeds 
while keeping design thrust constant, but also for the purpose of lowering 
unsteady propeller thrust levels. For a given propeller RPM use of more 
blades yields a higher value of reduced frequency associated with the blade- 
passage fundamental, and so in lower levels of noncompact blade-sectional 
loadings contributing to total thrust. 

2. An obvious suggestion is that the number of propeller blades be 
kept even in the pylon-mounted, once-per-revolution excitation mode. For 
this configuration odd blades would yield better coupling between odd 
harmonics of thrust and those of the inflow nonuniformity and therefore would 

be undesirable. By the same reasoning, if a pusher propeller is ever considered 
for general wing installation, it should contain an odd number of blades since 
the excitation would then occur twice per revolution rather than once. 

Wing Excitation by Wake of a Tractor Propeller 

1. Similar to Item 1 above, a large number of blades will yield higher 
reduced frequencies based on wing semichord and flight speed, and so, in lower 
levels of wing pressure distributions for this mechanism. 

2. Increasing the radius of the propeller should result in some decrease 
in strength of blade-tip vortices driving the downstream wing, essentially 
because keeping total performance thrust constant would then imply a more 
gradual spanwise load distribution for individual blades and a lower near-tip 
value of steady blade-sectional lift. 
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NOMENCLATURE 


b 

c 

o 

c 

w 

E ,E* 

g 

G 

h 

h 


blade semichord (Chapt. II) , wing semichord 
(Chapt. Ill) ? dimensional 

sound speed; dimensional 

wing chord ; dimensional 

Fresnel Integral and conjugate, Eqs. 11.25,26 

wing thickness distribution (g* nondimensional) 

function defined by Eqs. 11.30, 31 

wing camber distribution (h 1 nondimensional) 

instantaneous normal- to-wing vortex position ? 
dimensional 


k (R) 
m 

k 


k 

ac 

K 


l 


L 

max 

M(R) ,M 


N 

^mN ,P mN 


r 


R 


gust wavenumber at radius R; dimensional 
o)b/U, reduced frequency 

kcosA , normal-to-blade leading edge wavenumber, 
equal to k in part of Appendix A 

normalized acoustic wavenumber wb/c 

o 

Possio kernel 

arc length along blade midchord; dimensional 

maximum sectional steady blade lift ; dimensional 

circumferential Mach number for blade section at 
R, in Chapter II ? flight Mach number in Chapter III 

flight Mach number in Chapter II 

number of blades 

wing pressure distribution and spanwise transform; 
dimensional 

two-dimensional range measured from blade leading edge, 
Appendix A; nondimensional 

radial distance to blade section in question; dimensional 
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R* 


R t 


VV R eff 

R 


s 

t 

T 


m 

tot 


p tot 

"mN 


mN 





u 

m 

U 


w 


W / w 
mN mN 

i s 
w ,w 


x 


*b 


Radial distance defined in Figs. 11.4,5; dimensional 
blade-tip value of R* 

Tip, hub, and effective- tip radii in Eq. 11.28? dimensional 
R t /b 

local spanwise blade coordinate? nondimensional 
time? dimensional 

unsteady blade thrust, Chapter II? steady-performance propeller 
thrust, Chapter III 

harmonic of blade thrust; dimensional 

propeller thrust; dimensional 

harmonic of propeller thrust? dimensional 

harmonic of sectional blade thrust? dimensional 

nonuniform propeller inflow (Chapter II) , or tangential 
velocity for shed wake (Chapter III) 

normal-to-wing leading edge interference horizontal flows 
due to thickness and lift, respectively; dimensional 

normal-to-propeller disk components of u*,u*; dimensional 

m^ Fourier component of u; dimensional 

flR in Chapter II, flight speed elsewhere? dimensional 

flight speed in Chapter II ; dimensional 

propeller-induced downwash velocity; dimensional 

harmonic of w, dimensional? and its spanwise Fourier transform 

downwash velocity due to induced and shed-wake flows? dimensional 

normal-to- leading edge coordinate (Chapter II for blade, 

Chapter III for wing) ; nondimensional 

chordwise position of maximum camber displacement; dimensional 

blade-based coordinate, Fig. II. 3? nondimensional 

horizontal coordinate normal to wing leading edge, Fig. II. 2; 
dimensional 
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y(R*) #y t 
y 



x^ position of propeller hub 

blade mean-sweep function and tip value, Figs. 11.4,5; dimensional 
spanwise wing coordinate, Chapter III; nondimensional 
maximum camber displacement; dimensional 
yb 

x f + i z f ; dimensional 

normal- to-wing coordinate; Chapter III; nondimensional 
normal-to-wing coordinate; Chapter II; dimensional 
value of propeller hub 


a 

e 


wing angle of attack 

4-M 2 (R)cos 2 A (R) , Chapter II; or / 1-M 2 based on flight 
Mach number, Chapter III 
2 2 

/ 1-M cos A (Chapter II) 
f w 


6 


t 


Y 


angle of blade twist plus pitch 

strength of wake vortex sheet; dimensional 


r 

e 

A ( R ) 

w 

P ' 


P 


vortex strength; dimensional 
2 

p+k /6 in Appendix A; nondimensional 

angular position in blade polar reference frame (Appendix A) ; or 
circumferential angle in propeller disk (Chapter II) 
circumferential component of blade sweep, wing sweep 

wavenumber parameter, Eq. 11.24? nondimensional 
2 

equal to kM/B in Chapter III 

frequency parameter; Eq. III. 23; nondimensional 


background density; dimensional 
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wing thickness ratio 

perturbation velocity potentials? dimensional 
harmonic of $ 

wavenumber parameter, Eq. 11.23; nondimensional 
propeller rotational speed [rad/secl 
m AN 
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I. INTRODUCTION 

An important task of the present research effort is the prediction 
of some of the new wing or pylon forces likely to be introduced by the 
aerodynamics of high-speed propellers. Estimates similar to those reported 
here will serve as input to aircraft-specific models of structureborne 
vibration for the wing and fuselage until more accurate calculations become 
widely available. We consider two possible mechanisms for the generation 
of such airloads: 

First, chopping through the wake of an upstream supporting pylon 
will cause unsteady blade thrusts on a pusher propeller. If not properly 
isolated somewhere within the volume of the nacelle, the total thrust force 
will drive said pylon harmonically in a direction approximately colinear 
with that of flight, thereby producing fuselage vibration and cabin noise. 

The same can be expected of a wing-mounted tractor propeller cutting through 
the steady upstream interference field of the wing. The present study 
examines both modes of propeller installation. 

Secondly, under the assumption of zero external interference or cross 
flows, steady operation of a tractor propeller has an associated system of 
helical blade-tip vortices which appear frozen in the rotating frame of 
the propeller, but on the wing plane induce an unsteady downwash distribution 
that in turn gives rise to unsteady pressures and lift. 

Hanson^ - has developed what is probably the most general linear- regime 
formulation for lifting and nonlifting unsteady propfan flows. Given 
sufficient computing resources and time one could -apply his model to analyze 
rigorously both situations listed above (standard lifting-surface theory 
being a special case of his integral equation, which therefore could be used 
for the wing as well). As earlier stated, however, our objective here is 
the less ambitious one of temporarily front-ending our own fairly detailed 
structural models for wing and fuselage with rough estimates for the aero- 
dynamic forces which might drive them. To this end, predictions of unsteady 
thrusts for blades passing through either a pylon's wake or the wing's 
upstream interference field are here based on the following simple model: 
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The nonuniform propeller inflow field is broken up into a set of 

circumferential gusts at a chosen number of radial stations between 

hub and tip. For each of these we compute the value of blade sectional 

lift (thrust) using the approximation of local aerodynamic reaction 

2 

implicit in high-frequency strip theory. The model assumes zero blade 

pitch and twist, and so really addresses only the circumferential component 
of the aerodynamic problem ; it becomes more and more naturally valid for the 
radially outward, highly loaded blade sections at lower pitch. 

At every blade section we neglect coupling between steady-performance and 
unsteady loadings. However, Myers and Kerschen^ have determined that the omitted 
interaction could lead to 0(1) underestimates of airfoil radiated noise, 
so that their results rationally suggest that our predictions here could be 
low by 6 dB. On the other hand, since at each blade section we apply 
a Mach number based on only circumferential motion, the effective decrease 
in compressibility through neglect of the flight speed component should 
result in some overestimate of sectional thrusts, thus partially offsetting 
the above neglected effect. The model also starts out under the assumption 
of no interference among blades, but later demonstrates that this simpli- 
fication is indeed largely justified for the parameters in question. 

4 

The basic unsteady- aerodynamic theory is that by Adamczyk for a swept , rigid 
flat-plate airfoil passing through a gust? we apply it at each blade section 
using the local value of blade sweep, chord, and circumferential wave- and 
Mach numbers. Tip effects are simulated in an ad hoc fashion by allowing 
for decay of load amplitudes near the blade tip while keeping phases at 
their two-dimensional strip-theory values. In the final step individually 
integrated blade thrusts combine to yield the total thrust harmonic acting 
along the propeller axis. 

For the mechanism of wing excitation by the propeller wake, the model 
uses a system of straight potential vortices appearing on the wing's Trefftz 
plane as a set of points rotating about the propeller axis; the 
downwash induced is constant across the wing chord because the vortices 
are taken to be aligned with the freestream direction. The wing is again 
a flat plate of constant chord and infinite span. The unsteady-aerodynamic 
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theory this time is Amiet's^ generalization of that in Ref. 4 to include 
plunging motion, valid for high subsonic flight speeds and/or high 
propeller blade-passage frequencies (it should be pointed out that Amiet 
has developed his theory for chordwise downwash distributions other than 
constant, thus allowing for the possibility of modelling some of the effects 
of actual vortex tilting relative to the wing chord. This further 
generalization, however, is not pursued here) . 

6 

Weir has recently analyzed the low-speed case of this mechanism. 

Because the application to structurebome sound analysis for this 
mechanism uses continuous wing plating to interface with the source, 
predictions will here be in the form of a pressure distribution over the 
wing surface instead of a single load or moment, or running lift. 

In summary, Section II will be devoted to the prediction of the time- 
varying point thrust propellers could feel due to several sources of non- 
uniformity in inflow, while Section III analyzes unsteady distributed wing 
airloads induced by the steady vortex wake of a tractor propeller. 
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II. predictions of harmonic axial forces for pusher and tractor 

PROPELLERS IN NONUNIFORM STEADY INFLOW 

A. Introduction 

There are two prototypical configurations being considered for the 
new propfan technology. In one the propeller is mounted in the conventional 
sense, ahead of the wing, with the nacelle growing out of the latter. This 
mode of installation is currently undergoing tests on a Gulfstream II aircraft. 

In the other basic design each nacelle is held near the stern by an 
outboard pylon connected directly to the fuselage. This type of installation 
will replace the aft-mounted turbojets now on T-tail aircraft such as the 
DC-9. It uses pusher instead of tractor propellers, whose wake flows 
do not impinge directly on any part of the aircraft, which therefore should 
be free of that type of excitation. However, since each propeller must then 
cut through the wake field of the upstream pylon supporting it, new unsteady 
blade loads are generated and transmitted to the wing via the engine structure 
unless isolated effectively. 

This harmonic propeller thrust force is also present, though to a much 
less expected degree, for a tractor propeller passing through the upstream 
interference field of the wing. In the following sections we develop an 
analytical/numerical model to predict its spectra for both types of 
installation. 

A basic assumption in the development is that blade-blade unsteady 
interference is negligible. In Appendix A we use the blade loads so 
computed and verify that mutually induced flows and pressures are in fact 
an order of magnitude lower than corresponding quantities for the incident 
gust field. 


- 11 - 



B. Models for Incident Velocity Fields 


1. Pylon Wake Affecting Pusher Propeller 

The wake of a pylon will be modelled using available measure- 
ments for the wing of a typical large aircraft. 7 Fig. II. 1 shows mean-wake 
data for the 32%-chord station downstream from a wing's trailing edge. We 
shall construct the wake of a nonliftihg (symmetric) pylon by taking the 
average of velocity deficits above and below the wing plane, thereby 
generating a figure similar to II. 1 but for 

„ [u(z + ) +u(z ) ] 

1 2 — f * 

We shall also adjust the data to generate wake profiles at stations other 

8 

than the 32%-c^ downstream point. The procedure will be to apply Schlitchting' s 
simple mixing- length model for the diffusion of the two-dimensional wake 
behind a cylindrical body. He calculates that the rates of widening for the 
far wake and of reduction for the maximum velocity deficit are respectively 
proportional to the square- and inverse square- root of distance downstream 
of the object. Thus we have that u(x,z=0) /u(x=32%c^,z=0) can be obtained 
from / 1. 32c w /x and b^ (x) / b w 32%- c from ^x/1. 32c w , where b^ here denotes the "wake 
width" and x the distance downstream measured from the leading edge. By 
referencing x to the leading edge we generate a wake profile still fairly 
sharp at, say, one full chord downstream of the trailing edge (x=2c^) . If 
we used instead the trailing edge or even the mid-chord point, a more diffuse 
profile would result. The sharper profile is chosen because it yields 
conservative estimates for the blade-passage and higher harmonics of total 
propeller thrust. 

2. Upstream Interference by the Wing 

In the case of a tractor propeller the upstream wing inter- 
ference for lifting and non- lifting wing flows at cruise will be computed 
by applying potential theory to the actual airfoil section of interest. The 
model uses a wing of infinite span and thus ignores wing-tip and fuselage 
effects at the propeller location. Compressibility and wing sweep enter 
the calculation through the Prandtl-Glauert factor /l-M^cos 2 ^ w ' M f cos ^ w be i n 9 
the component of flight Mach number normal to the wing's leading edge. 
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Fig. II. 2 shows the idealized wing propeller system for the upstream 
influence calculation. All complicating structures have been removed for 
simplicity. The perturbation velocities u* in the direction normal to 
the wing's leading edge are readily found from classical potential theory, 
e.g., Ref. 8: to be, for thickness and lifting flows, 


★ 


u 

t 


U _cosA 
f w 




(x f -x)g' (x) 
—— — - 

(x f -x) “+Z 


(II. 1) 



where c^ is the wing chord, the flight speed, A^ the wing sweep, and 

Z =z_3_, where 0 = / l-M^cos^A , with M =l! fi /c , the flight Mach number, 
rfr f f w ffo 

In Eq. II. 2 Y=x f +iZ f . The f subscripts on x f ,z f , and denote field points; 

on f denotes the flight speed of the aircraft. In Eq. II. 1 T is the 

airfoil's thickness ratio and g' (x) is the surface slope function from 

leading to trailing edge. We shall model our wing as a four-digit NACA 

9 

airfoil, for which the thickness function g* is given, by 


/C 

g 1 (x) = .7423 / -.6300 - "3.5160 (~) 

l "w 

2 3 

+ 4.2645 (— ) -2.030(— ) (TI.3) 

c c 

w w 

For the above "real" airfoil Eq. TI..1 is inapplicable if one's objectives 

★ 

are pressures (or u ) in the vicinity of the leading edge, essentially 
because the round nose violates the smal 1 -disturbance assumption of the 
theory. This defect, however, is local and should not seriously affect the 
interference calculation of interest here. Loading solutions for an 
elliptic airfoil discussed by Van Dyke 10 , for example, indicate that even 
though the small-disturbance solution fails at the leading edge, it 
nevertheless yields the correct interference several nose radii av/ay . 
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In Eq. II. 2 h' is the slope of the mean-camber line, for a four- 
digit airfoil given by 


2y (x -x) 0 < x < x 

a a — — a 


h' (x) = 


"2y (x-x ) x = < x 1. c w 
a a a — — w 


(c -x ) 
w a 


(II. 4) 


where x is position between leading and trailing edge where the camber 
a * 

displacement has its maximum value of y^. In Eq. II. 2 a stands 
for angle of attack. 


As we have said, Eqs. II. 1 and II. 2 give the perturbation velocities 

in a direction normal to the swept leading edge. Of interest here are 

corresponding quantities u , u in the direction normal to the propeller 

* t x, 

plane; in terms of u these are obtained from 
* 

u = u cosA (II. 5) 

t t w 


★ 

u rt = u rt cosA (II. 6) 

i i w 


Referring to Fig. II. 2, we note that field points on the propell er disk 
and z^ are given by 


x* = x +R sinA cos0 , x <0 (II. 7) 

f Of w o f 

Z-. = z +R sin0 (II. 8) 

f °f 

where the value of R ranges from R^ at the propeller hub to at the 
propeller tip. 

Substituting Eqs. 11.7,8 with 11.1,2 into Eqs. 11.5,6 yields the total 
upstream interference field 


The percent chord values of y ,x take up the first two digits of the 

4-digit NACA designation, and ihe a thickness ratio the last two digits. Thus 

a NACA 2410 airfoil has y /c =.02, x /c =.04, T=.l. 

aw aw 
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(II. 9) 


u(R,e) = u t (R,e)+u^(R,e) 

normal to the propeller plane as experienced by the blade section at radial 
position R as it travels from 9=0 to 6=2 tt. 

C. Unsteady-Aerodynamic Model for Propeller 

1. General Coupling Between Velocity and Lift 

The nonuniform inflow field at the propeller plane, whether 
caused by upstream interference or by the pylon's wake, is now expressed at 
every radial station as a superposition of circumferential gusts: 

oo 

u(R,0) = u (R)e~ im0 (II. 10) 

m=-°° m 

2ir 

u (R) = f ^u(0)e lmS (II. 11) 

m J 2tt 

0 

On the rotating frame of the propeller the speed of travel for each 

of the above gusts is U ( R) R , R <j*^_R ; Q is blade rotational speed in radians/ 

n— 

sec, and the phase factors in Eq. II. 10 become in(0.-fit) for the j propeller 
blade turning in the clockwise direction (9 is defined positive in the 
counter-clockwise direction) . Interpreting angular displacements in terms 
of local arc length, we have that 

m0 = ~ # R0 = k^(R) b [cosA (R) x^+sinA (R) • s] (11.12) 

where k (R)=m/R becomes the gust wavenumber value at radius R , A ( R) is the 
m 

local blade sweep angle, x^cosA^ the blade-based origin of the locally 
rectilinear circumferential direction , and s is the local spanwise, or radial, 
coordinate. The situation for a typical radial station is portrayed in Fig. II. 3. 

If the blade sweep's "radius of curvature" is much greater than the 
acoustic wavelength everywhere between hub and tip it should be possible to 
analyze each radial blade position as a locally straight planform passing 
through a gust of wavelength 2fTR/m convected at local Mach number M(R)=QR/c o 
at local blade sweep A(R). Fig. II. 4 shows how the latter is geometrically 
defined and calculated. We shall consider a family of radial sweep distri- 
butions given by the relation 
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(11.13) 


y(R*) = y fc (R*/R*) 


where 


R = v^* 2 +[y (R*) ] 2 , R*= / R^-y^ (II. 14) 

where as indicated in Fig, II. 5 y^ is the in-plane position of the blade 

tip relative to that were the blade not swept; y ^=0 in the analytical 

model corresponds to unswept blades. From Eq. 11.13 and Fig. II. 5 we have that 

A(R)=A(/R* 2 +y 2 ) = tan" 1 -* ) ] -tan" l fy/R] (11.15) 

3 R* 

We temporarily postpone further discussion of details of the unsteady 

aerodynamic calculation and for now simply observe that its goal is to 

yield a value of modal thrust T for each of the N blades cutting through 
th ^ 

gust m. For the j blade then, the time history of unsteady loading would 
be given by 

E -im(0 . -Q t) 

T e 3 (11.16) 

m 

nF - 00 

where T is proportional to u . Taking 0 . =2tt ( j-1) /N , j=l,2...,N for N 
m m j 

equally spaced blades, we obtain the total propeller thrust as 
tot N 

T(t) = £ T. (t) (IT. 17) 

j-1 3 

Substituting Eq. 11.16, interchanging orders of summation, and noting that 


N -im2iT ( j-l)/N 

£ e 

j = l 


0 if m^multiple of N 
N if m=multiple of N 


we have that 


00 

T(U = £ T tot 


T t0t =NT 


(11.18) 


(11.19) 



obtained 


where Nft is the blade-passage frequency; T^ is proportional to u^, 
from Eq. II. 11: 


2tt 


",* (R> 


r d9 

•J 57 U(R - 


0)e 


imN0 


(11.20) 


Quantities NT „ for m=l-,2..., represent the propeller thrust harmonics 

of primary issue here. The corresponding gust wavenumbers of interest are 

k (R) =mN/R , etc. 
mN 

2. Calculation of Modal Blade Sectional Thrusts 


The sectional thrust (Ib/ft) at radial station R corresponding 

to gust wavenumber k kf is given by (gust interaction- and sweep angles both equal A) 

mN 


27rbpUcosAu 


■ik s 
s 


T m iR) r 


mN 


/ k b cosA/0 +p 

mN 


l^i E*[2(y-i|>)] + £(tl ^) for M>t 


1-i EJ2QMO } 

* /- t 

/ if#— H 


+ f ( y , i|y) for )K i|j 


71 1 71 


- 1 - + - (1+ ~~) E* ( 4(i ) 
ir 8p 


. -i4u 1 

- i® 

3/2 yr-J 

7T v'2|l J 


for 


(11.21) 


with 




i 1_ T -2i (p-i(<) ] /TiT 

3/2 y-i|y L® \* I* H l> 


2 li E* [2 (y+i|i) ]- — } 
||H|> ^ ) 


+ - E*(4y) J (11.22)* 

* The last two terms of Eq. 11.22, -E*(4y), appear misprinted in Ref. 5 

as i|i-E*(4y). 
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and where 


* = 


2 3 

k bM cos A 
mN 

e 2 


(= k M 2 cos 2 A/6 2 ,k =k b cosA) (11.23) 
x x mN 


P = 


/ „2 2 * 2 A 

/M cos A-sin A 

e 2 


(11.24) 


with A=A (R) given by Eq. 11.15; b=b(R) is the radially varying blade 

2 2 

semichord; again, k =mN/R, M=M ( R) -U ( R) /c =QR/c and 3 = / 1-M cos A; p is 

mN o o 

background density. E ( a) denotes a Fresnel integral of argument a; 


E (a) 


-/ 


dt e 


it 


(- 


1 + i 


as a-><“) 


(11.25) 


0 / 2irt 

so that its conjugate E* is given by 


E*(a) = 


d 

-f 


dt e 


-it 


(- 


1-i 


as a-**) 


(11.26) 


/ 2fft 


-ik s 

s 

The exponential factor e in Eq. II. 21 keeps track of the relative phase 

of sectional loads for different positions along the blade radius. Referring 
to Fig. II. 3 for the straight blade in rectilinear flight, the quantity 
s . sinA denotes the flow direction projection of the spanwise distance 
between two blade points. The phase factor should therefore be f k^b. s . sinA} 


or exp( -ik s} with k =k bsinA 
s s mN 


In the standard literature the spanwise 

coordinate s usually carries the symbol "y" so that exp{-ik^s} is 

expf -i M k y"} , with "k y"=k bsinAy in agreement with the above argument, 

y y mN 

The exp { -ik s} factor will be simulated here by analogy to the 
s 

rectilinear case. The arc length difference in the direction of the in- 
coming circumferential gust between radial station R and that at the hub is 
Rsin \/R, with y given by Eq. 11.13. The phase factor therefore becomes 


-ik s 
s 


-imNsin y/R 


(11.27) 


The slight relative shift in center of pressure, or lift (1/4 chord for 
steady lift), due to the variation in circumferential Mach number and gust 
wavenumber from blade hub to tip will not be taken into account (cf. discussion 
in Ref. 4 , p. 24) . 
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Finally, we wish to modify the basic sectional thrust distribution 


to account for the / R^-R load dropoff as R approaches the tip. Thus 

we choose a value of R=R where the sectional thrust will achieve 

err 

approximately its maximum value, say at 

R eff * \ + 7541 VV <1I ' 28 > 


and define the total blade thrust harmonic T M for use in Eqs. 11.16,19 

mN 

as follows: 


mN 


■I 


dl T [I (R) ]G(R) 
mN 


(11.29) 




where l denotes length along the blade and 


G(R) 


for R, < R < R 

h — — eff 


.(a +a R) / R -R for R .. < R < R 
\ o 1 t eff — — t 


(11.30) 


(11.31) 


with 


a 

o 


eff 


2 ( V R eff> 


3/2 


/R -R 
t eff 


(II. 32) 


a l = 


2 ( V R eff ) 


3/2 


(11.33) 


Constants a ,a, are determined from the continuity requirements for the 
o 1 

sectional load and its slope at R ... 

eff 
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D. Application of Model 


1. Propeller/Aircraft Geometry 

Figures II.6a,b show front views respectively of a Hamilton 

* 

Standard propfan prototype and of one of its blades. The indicated di- 
mensions are approximate projections on the propeller plane of actual 
propeller chord widths, etc. Our calculations will be based on them. 

The blade sweep distribution from hub to tip uses the midchord curve 
shown in Fig. II. 6b; in applying the model earlier described and depicted 
in Fig. 11.4,5 we found a curve. with the power n=3 adequate for the 
present case. 

Fig. II. 7 shows top and side views of the installation geometry for 
a Gulfstream-II aircraft, and dimensions drawn for our calculations. 
Unfortunately, we were unable to obtain similar schematics for a pylon- 
mounted propeller; our calculations for this case are here based on 
assumed propeller locations downstream in terms of pylon chords. We use 
the 32%c^ and 100% c^ wake positions as measured from the pylon's trailing 
edge. 

2 . Models of Incident Velocity Fields 

(a) Wake Signatures for Wing or Pylon 

Fig. II. 8 shows the result of fitting a two sided 

Gaussian to the wake measurements shown in Fig. II. 1 for a lifting 

aircraft wing, and modelled velocity deficits experienced by two radial 

points on the propeller in their circumferential trajectory. The figure 

applies to a propeller with axis lying on the plane of the wing, which 

extends in the spanwise sense on both sides of the propeller. Only the 

0=7T spike will be used for a nonlifting pylon, and then symmetrized 

about its maximum value. The signature felt by the blade section near 

the hub (R, ) is wider than that at R --=.8 R . because the effective wake 
n eff tip 

takes up then a greater fraction of the total trajectory. 

As a further illustration, Fig. II. 9 shows the situation for a pusher 
wing-mounted propeller with axis at a distance above the wing plane 
corresponding to that of the Gulfstream-II nacelle (although a pusher 

* I. Loeffler, NASA Lewis Research Center, personal communication. 
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wing-mounted propeller is presently not a candidate concept, we examine 
it also for completeness) . 

(b) Upstream Interference for Wing 

Fig. II. 10 shows predicted interference wing flows for 

the circumferential trajectory of R=R _ ==.79 R . -3.5 ft. The top insert 

eff tip * 

shows both the cross-section of the^ Gulf stream-II wing in a plane parallel 

to the flight direction, arid the equivalent 4-digit NACA airfoil section in 

a plane roughly perpendicular to the leading edge (note Fig. II. 7)? 

y /c -.013, x /c =.19, thickness ratio t= 10%. A value of 1.5° was taken for 
aw aw 

the angle of attack. Point 9=0 is inboard and axially closest to the wing 
(middle sketch and Fig. II. 2a). Nonlifting interference is always negative? 
it is maximum at 9=0. Lifting interference is positive above the wing 
(because the lift is positive) and negative below; it vanishes identically 
at the two points where the cuts the extended plane of the wing. 

3* Frequency/Mach Number Regime of Validity 

Fig. II. 11 shows the sweep distribution of the midchord 
curve traced out in Fig. II. 6b. The other curve is a plot of Eq. 11.24. 
After comparing approximate results v/ith other more rigorous calculations, 
Amiet (Ref. 5, p. 1078) has determined that Adamczyk's high-frequency 
theory for a leading-edge flat plate passing through a gust at zero angle 
of attack should be reasonably accurate for u>ir/4. If trai. ling-edge effects 
are included (as we have done here) , he states that the results apply 
beyond a much lower value of vj . As far as compressibility and frequency 
effects in the propeller plane are concerned, we conclude that the \i curve 
in Fig. II. 11 confirms the need and relevance of the noncompact aerodynamic 
model being applied. The values shown correspond to the fundamental 
N 0=226 Hz (m=l , N=8) . 

4. Predicted Blade Loadings and Total Propeller Thrust 
(a) Pusher Propellers 

Figs. II-12a,b are for a pusher propeller on an inboard 
nonlifting pylon at 32% and 100% chords downstream of the trailing edge. 

The axis of the propeller coincides with the midplane of the pylon. Each 
figure shows two curves: the wake field's fundamental harmonic as a 


* W.E. Arndt, R.W. Nazarowski , Lockheed Georgia Co., personal communication. 
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function of radial position, u ( R) =u (R) ; and the corresponding sectional 

mN N 

blade thrust as given by the integrand of Eq. 11.29, The 32%-c^ downstream 
wake is only slightly richer in Fourier component than the 100% -c^ wake, 
with the largest differences appearing at the hub station. 

tot 

Figs. 11.13, 14 show predicted harmonics of total thrust 

(Eq. 11.19) for the two cases of Fig. II. 12a, h. The 32-dB level for the 

fundamental corresponds to a 40-lb force acting along the propeller axis. 

Major differences between 32%-c and 100%-c curves appear to be postponed 

w w 

to higher frequencies. 

Figs. II. 15a, b show results for the unpursued concept of a wing- 
mounted pusher propeller; they correspond to propellers with (a) axis on 
the wing plane and, (b) axis elevated 1.78 ft above the wing plane. The 
axial propeller location is 32%-chord downstream from the trailing edge. 

The roughly doubled values in Fig. II.15a relative to those in Fig. II. 12a 
reflect mostly the more effective coupling that exists for the "two-per- 
rev H wing over the once-per-rev pylon, given that the number of propeller 
blades is even (8). Fig. II.15b for the raised propeller displays a fairly 
erratic behavior in fundamental harmonic of incident velocity field (and 
associated sectional thrust) . 

Figs. 11.16 and 17 show corresponding predictions of total thrust 

harmonic content. The 38-dB value for the fundamental in the z =0 case 

° f 

corresponds to a 70- lb force. The reduced value for 226 Hz in Fig. IT. 17b 
can be attributed to the relatively poor coupling some radial stations 
have to the 8 th Fourier modal number (Fig. II. 15b). 

( b ) Tractor Propeller on a Gulfstream-II Wing 

Fig. 11.18 shows the propeller airload response to the 
combined lifting and nonlifting flows of Fig. IT. 10. The extremely low 
values (-48 dB=4xlO ^lb) suggest that this mechanism will not be significant. 
The reductions in level have two causes: first, the maximum interference 

amplitudes are low relative to those for the wakes previously examined, say 
2-3% vs 20-25%; this difference, however, would account for only 20 dB 
relative to the wake results. More significantly, of course, is that the 
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propeller is sufficiently removed in the upstream direction that the 
interference signature is a very smooth function of 0 at all radial 
stations. Since the response couples only to multiples of the funda- 
mental Fourier counter, which is high due to the number of blades, it 
would appear that low levels are virtually guaranteed for this mechanism. 

E. Predicted Benefit of a Large Number of Blades in Reducing Harmonic 
Propeller Thrust 

Eqs. 11.19 and 21 point out that the theoretical total unsteady 

propeller thrust is proportional to p(ftR)w (k b) ^ .N . Blade Area, where 

mN mN 

k E=mNQb --/(ftR - J =mNb . _/R an effective or average reduced frequency. 
mN err eff eff erf 

We have assumed that its value is sufficiently high to approximate Eq. 11.21 

by its high-frequency asymptote, i.e. , with E*=(l-i)/2, etc. If we now 
eff eff 

write w =a fiR, where a is an effective instantaneous anale of attack 
mN mN mN 

in the quasi-steady aerodynamic sense, it follows that the unsteady thrust 

2 1 

harmonic for the propeller is proportional to p ( ftR) . N. Blade Area. 

Now, the corresponding quantity for the design steady thrust may be 

2 2 2 - - 

similarly written as p [U + Q R ] .N. Blade Area.c , where c is an effective or 

f L L 

spatially- averaged lift coefficient for each blade. With QR/U f =0(l), 
we therefore obtain that the ratio of harmonic propeller thrust to propeller 
design thrust is proportional to 


T 

T 


tot 

mN 

tot 

design 


1 


(k J j)c r 
mN L 


1 

N 


(IT. 34a) 


(TT . 34b) 


a vanishing quantity for high N for a family of propeller designs which keeps 
performance thrust constant. Having a large number of blades therefore makes 
sense from the perspective of propeller unsteady aerodynamics. 
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III. PREDICTIONS OF WING AIRLOADS INDUCED BY A PROPELLER WAKE 


A. Introduction 

The model for the propeller wake will here be essentially the 

same we developed in Ref. 11. A system of rotating straight potential 

vortices , one for each blade, radially located at the propeller tip 

positions. The path, or position, of the vortices is assumed unaffected 

by the presence of the wing — unlike in more sophisticated treatments 

12 

by Hardin and Mason, of a similar situation at low speed, and by Lindblad, 
for high speed. Also, we take each vortex to be aligned with the freestream 
direction, so that induced velocities over the wing chord are constant. 

Even for this simple representation the system will be seen to induce 
a fairly complicated unsteady downwash distribution as function of spanwise 
location. In Ref. 11 we circumvented the development of a theory that 
would couple wing sectional airloads due to such an inflow by arguing that 
at high speeds (M=.8) and reduced frequencies individual spanwise stations 
should react locally, in a strip-theory or lifting-line sense, because 
noncompactness effects should then tend to concentrate unsteady airloads 
at the wing's leading edge. In the present study we remove this approxi- 
mation. Now different wing sections are allowed to influence each other, 
and the computed airload at each point on the wing surface depends on that 
at every other point. The wing remains unswept and of constant chord. The 
model is applied to the Hamilton Standard propeller/Gulfstream II wing 
combination. 

B. Propeller Wake Model and Wing Aerodynamics 
1. Analysis of Downwash 

Figure III.l shows the wake model. Each vortex contributes a 
downwash field on the wing plane given by 


2tt 


Y-Y 


(y-y Q ) +h 


2 

o 


with y =R u cosfit, h =R sinftt; so that the total influence of the N straight, 
o t o t 

equally spaced vortices is 
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w(t,y) 


(TII.l) 


£ ^ y/R-cos[^t+2TT(k-l)/Nl 

OttD 1 -2 - 

t l+(y/R) -2(y/R)cosrnt+2ir(k-l)/Nl 


where y and R stand for ?/b,R t /b, with wing semichord b assumed constant. 

The analysis to obtain the corresponding harmonics begins with an 

i 8 ia . 

identity involving two arbitrary complex numbers x=e , z=re (e.g. , 

Ref. 14) . 


-L- + — 

T-Z * * 


1 -r 


2 


(III. 2) 


t -z l+r“- 2 rcos (a- 0 ) 
where * denotes conjugate. 

For r<l we may now express the left hand side of Eq. III. 2 as 
00 

1 O V fsv-Q\ 

x * x _ ■ l - m y 

m=l 


Since r, the magnitude of z, is by definition positive, this relationship 

is immediately applicable to the range CKy/R<l. Thus, letting y/R=r, Gt=a, 

-3 = ' and performing some algebra one finds that 

N 


y/R-cosfftt + 2lT ^ ^ ] 


1+ (y/R) 2 -2 (y/R) Cos [Qt^ 


„ , ,- % m-l 2 iT(k-l) 

= - Z (y/R) cosmf^tt 


(TIT. 3) 


By a similar calculation for r>l the left hand side of Eq. III. 2 takes 
the form 


. I 


« cosm(a- 8 ) 


\ 1+2 E. 
J m=l 


m 


and for l<y/R<°° we find that 


_ 2tt (k-1) , 

y/R-cos[ftt + — ] 

N 

- 2 3 2tt fk-l) 

l+(y/R) -2 (y/R) cos [ Ot + 1 


00 cosmtat + 
E rt — 

m =0 


2 n(k-l) 

N 


, y - x m+l 
(y/R) 


(III. 4) 
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If the cosines on the right sides of Eqs. III. 2 and TTI.4 are now expanded 

out, and the sum E is interchanged with E, the following two factors appear 
k m 


. £, cos [ 
k=l 


27Tm(k-l) 

N 


k£i sin[ 


27Tm(k-l) 

n 


the first of which is zero unless m is a multiple of N , and then it equals N. 
The second sum is identically zero. These conclusions may be used to find, 
for example, that the right side of Eq. III. 3 becomes 

-N £ (y/R) Nm ^cosmNQt 
m= i 


The simple manipulations used to obtain the series expression of w for 
0<y<°° may be adapted to the negative ranges -l<y/R<0, -°°<y/R<-l. One writes 


1-r 


1-r 


2 2 

1+r -2rcos(a-f3) 1+r -2 (-r) cos (a-0+iT ) 


„ 2 tt (k— 1 ) , 2 tt (k— 1 ) . . ^ 

and then replaces y by -y>0 and - by turn the previous 

results . 

The final series representations for w=w ( t ,-°°<y < °4 are found to be 


w(t,y)= 


TN 


4'irR 


» imN Qt 
r e u 

£ m 

m=-« 




- N|m|+1 

) 


R<y <c 


i i , “ M I mkl 
■ , , . N m +1 . R . 1 1 - 

V (-1) 1 1 ( ) -co<y<-R 

~Y 


and 


00 imNftt , ~ . 

~ e (5 -1) 

£ om 

m--°° 


(y/R) N ' m 1 0<y <R 


• 1 ) N l m l -1 ( - y/ R) N l m l 1 -R<y <0 (III. 5) 


where v =2 or 1 for m^O, m^O, respectively; Kronecker delta 5 equals 1 
m om 

when m=0 and is zero otherwise. Identifying w T (y) in the above, we have 

mN 

that 
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/ V TN 
W mN Y 4 ttR 


m 


(R/y) N ^ m ' +1 R<y < < 


, N | m | + 1 . R . N|m > +1 - 

(~1) I I (— ) — co<y <-R 


(5 - 1 ) 

mo 


(y/R) N m 0<y <R 


, _ . Nm-1 . -Nm-1 - 

(-1) 1 1 (-y/R) 1 1 -R<y <0 


(III. 6) 


As discussed in Ref. 11# it is of interest to note that the harmonic 
decomposition in Eq. III. 5 yields a nonzero value for w Q (|y|>R). This 
11 static" component of downwash would affect the steady performance of the 
wing. 


We shall need the spanwise Fourier transform 




00 

/ dy elXVw mN 


(y) 


(III. 7) 


of the above ; we obtain 


^mN U) = 2 ^ exp{ ¥ < N M +1)} 


[ 


V R N l m l +1 

m 


(N m ) ! 


! N ] m | (N j m| — j ) !cos [XR- ^■(N|m|-j + 

j-1 ^ 


-j + 2) 1 


-jYj- si ( | X | R) .X N ' m ' 


<1~0 (N|m|-1) ! (N|m|-1 ( -N | m| - j f _ n j + 1 


mo 

-N I ml -1 
R 


( N|m| 

i & 


— T — cosfXR- ^(N|m|+j-l) 

X 3 { N I m I - j ) 1 ’ 


sinXR ! 

+ . 


(III. 8) 


15 


where si(z) is defined as 
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si (z) 


00 

- 7 dt 


The quantity in large square brackets is real; it is odd in X for even 
values of mN and even in X for odd mN. Since the exponential factor outside 
equals (-1) (N M +L)/2 if Nm i s odd and i (-1) N I m l for even Nm, it follows 
that in general - . 


w = w k ,< x ) 

mN mN 


(III. 9) 


where * denotes complex conjugate, in agreement with the observation based 

on Fig. III. 2, which shows 4nR w (y)/TN (i.e. , m=l) for N even or odd, that 

t N 


i /dy sin(Xy)w (y) for Nm even or 
/ mN 


(III. 10a) 


w ( X ) = 
mN 


dy cos(Xy)w m (y) for Nm odd 
* mN 


[III. 10b) 


i.e., that Eq. III. 10a is odd in X while III. 10b is even in X. 

2 . Aerodynamic Response 

The harmonics of the downwash distribution of a single 
propeller have been shown to be always odd about the propeller axis if 
the number of blades is even, and either even or odd (depending on whether 
Nm is odd or even) when the number of blades is odd. 

The boundary-value problem for the perturbation velocity 
potential <t>(x,y,z,t) on the frame of reference of the wing is 


l jL + jL + jL\ , b 2 d 2 $ _ 

\ 3x 2 3y 2 3z 2 / c 2 Dt 2 


0 for all x,y,z 


(III. 11) 


— (0<x<2,y,o,t)=-w(t,y) 


(III. 12) 


■( x>2 ,y ,o,t) =0 


( 1 1 1 . 1 3 ) 


where c is sound speed and D/Dt=3/3t+(U/ fc) 3/3x;x,y , z have normalized by wing 
o 

semi chord b. We now let 
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oo 


«x,y,z,t) = 1 ] KI (x,y, z)e lmNQt 


m=-oo 


mN 


(III. 14) 


as required by Eq. III. 5; also, with u)=mNft # k=wb/U , and p=kM/S , we 
let 


. , . * , . i|iMx 


mN 


mN 


(III. 15) 


and obtain that 

2 * 2 * 2 * 

a $ a $ a <*> 

mN mN mN 2 ,* 

3x 2 3 (By) 2 3<ez ) 2 


(III. 16) 


- * , . -iyMx 

* mN (0<x<2 ,y ,0) __ W mN (y)e 


3(Sz) 


6 


(III. 17) 


Finally, letting Bz=Z and 


00 

*nV X ' y ’ Z) '/ 


dX e" lAy $* ( x , Z ; X ) 
— mN 

2 it 


(III. 18) 


to conform with Eq. III. 7, yields 


2 -* 2 -* 

3 4* 3 ‘1 

mN mN -2 7 * ^ 

— ~ + 7 — + h 4 KT =0 all x , Z 

3x 2 3Z 2 


( L r.1.20) 


~ r* w -iyMx 

% (0<x<2 ,0) ^ 8 

IT - B 


(III. 21) 


(^4 + f)i* (x> 2 ,°)=° 

a 2 3x mN — 

0 


(III. 22) 


where 


/ 2 .22 . 1 . 2.2 2 
y = / y ~X /B =-1 / X /3 -y 


(in. 23 ) 
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the latter equality actually being a consequence of the choice of branch 

2 

cuts for oj> 0 to satisfy the radiation condition (for p >0 is assumed to 

- 2 

have a small negative imaginary part which becomes large when y <0) . 

Since the propeller wake model uses straight vortices aligned with the 
flight direction, the induced downwash is constant across the wing chord. 

We now adapt an approximate solution by Amiet 5 for an airfoil in oscillatory 
plunging motion at high flight speed. Since Amiet has considered only the 
purely two-dimensional constant downwash problem , the e*xpression 
for loading given below constitutes a generalization analogous to Adamczyk’s^ 
skewed-wing analysis for a convected gust. We write the final result for 
the y-transform of the pressure distribution on the lower surface of the wing 


p % (0<x<2; z=0 ; X) » 
mN 


pUw (X) 

mN (1+i) 

/72 2 2 

/ y -y M 


{ kE* [ ( y-yM) xl - i / \ — expf-i (y-yM) x} 

[ v 2 irx 

+ E*[(fi4vM) (2-x) ]- (~) | (III. 24) 

from which we compute the surface pressure everywhere 


P mN (x ' Y '° 


QO 

-f 


f4 1>y P«. u ’ 0 ' 

2v 


; X) 


(III. 25) 


purN (l-t-j ) 

6R t 2, 2 


00 
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dX 


J-2 2 2 

/y-yM 


kE* [ (p-yM) x] -i / expf-i (p-yM) x} 
y ^ttx 


+ E*[(p+yM) (2-x) l-(~) 


• Re { exp (^-(N (ml +1) -iXy ) w < X ) } 
i mN 


(III. 26) 


- 30 - 



where w (X) stands for the quantity within large square brackets in Eq.III.8. 
mN 

Since p (-z)=-p „(+z) the actual loading, defined in the positive-lift 
mN mN 

sense as p (0-)-p (0+) , is therefore given by +2 times Eq. III. 26. Also, 

mN mN 

for m<0 p =p* , where asterisk denotes conjugate. 

mN (-mN) 

The large-y behavior may be found by deforming the contour impliait in 

Eq. III. 25, in X space, from the real axis to the several branch points present 

in the integrand; thus, for example, for y>0 one could deform to the lower half 

2 

plane to pick up contributions at X=+p0 =kM=k , the normalized acoustic 

2 22 - 1/2 aC 

wavenumber wb/c (from the (p -p M ) factor) and at X=+p0 (from a branch 

o 

point of p). It is obvious that of these the X=p0 point is clearly a weak 

contributor compared to X=k . As | y | -*» a straightforward asymptotic 

ac 

analysis therefore indicates that p ^exp{-ik lyl}, a net acoustic wave 

mN ac 1 1 

having its "origin" within the area of strong aerodynamic action. 

It is important to point out that solution III. 24 violates the condition 
of zero pressure on the plane of the wing both ahead of the leading edge and 
along the wing wake; it is the result of constructing a preliminary solution 
which satisfies the downwash boundary condition over the entire plane, and 
then correcting, by cancellations, pressures both ahead and behind the wing. 

As discussed in Ref. 5, for the chordwise-constant downwash case these 
corrections are found simultaneously rather than sequentially, with an out- 
come that does not satisfy pressure continuity either forward or aft of the 
wing. More corrections would result in better approximations to the computed 
loading. 

3. Estimation of Vortex Strength 

The estimate of strength for blade-tip vortices is here based 

on the following model: We assume that all blade-bound vorticity 

beyond a chosen radial station rolls up into a single concentrated vortex 

(Fig. III. 3). We may take this radial station as that where the sectional 

loading is maximum, although for our purposes even this is unnecessary. 

Assuming a simple triangular shape for the lift with a prescribed maximum 

value L , one obtains that 
max 
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(III. 27) 


2 L max R t V cosB t 

regardless of radial point where the maximum is assumed to occur; in the 
above T is total steady propeller thrust and 0 is an effective or average twist 
plus pitch angle. The term in parenthesis denotes each blade's contribution 
to total performance thrust. 

Applying these ideas at the tip radial position, we next invoke 
the lifting-line relationship, 

pr/u 2 +(ftRj 2 = L (III. 28) 

t max 

from which finally we obtain that the leading factor in Eq. III. 26 is 

P urN _ JL . 2 . U (III. 29) 

R t r 2 cos0 t /u 2 +(nRt ) 2 

suggesting that, for a given propeller design thrust T, lower pressures are 
obtained if the propeller radius is increased. 

C. Application to Gulfstream-II Wing 

For the Gulfstream-II we take b=7 ft (14-ft chord at engine location) 

o 12 

and the total propeller thrust T at cruise as 1460 lb; [3^=57 • Flight Mach 
number M= . 8 ,* density and sound speed values correspond to an altitude of 
35,000 ft. The results discussed below are for the fundamental frequency, 

226 Hz, which at 35,000 ft has a wavelength of 4.3 ft. 

Figures III. 4-7 show predicted pressure distributions p^(x,y,z=0 ) 
over chords at y=R/2, R, 1.5R, and 4R. All results display the l//x leading- 
edge singularity. For most spanwise stations the singular region appears to 
be quite confined in extent, following one's expectations given the high 
speed/frequency of the situation (p=28) . As earlier discussed, the nonzero 
values of loading computed along the trailing edge are an artifact of the 
approximate lifting-surface theory. 
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Figures III. 8-10 are spanwise plots of p (x,y>0 2 ^o~) for chord 

N f 

stations x=.5, 1, 1.5 (quarter, half, and three-quarter chords). For y<0 the 
results would be inverted mirror images of those shown, due to the antisymmetry 
of loading about y=0 for an even number of blades. Point 4R corresponds 
to a dimensional distance of 18* outboard of the propeller axis. The 
numerical solution displays the acoustic wavelength for large y in agreement 
with the discussion following Eq. nr. 26. The decay rate is l/|y|*. 

The most interesting feature of Figures III. 8-10 is the smoothness 
of the aerodynamic* response near y=R to an input downwash which is dis- 
continuous there (cf. Fig. III. 2a). The loading is generally nonzero at 
y=R because y=R is not a point of spanwise downwash antisymmetry. Neither 
of these observations would apply to the loading obtained from the strip- 
theory model, which requires that the spanwise behavior of the pressure be 
the same as that of the input downwash, which, incidentally, lacks the 
acoustic wavenumber at large distances from the origin. The input downwash 
is strictly hydrodynamic in nature, while the aerodynamic response to it 
should be definitely wavelike, or acoustic. 

We have noted that for an even number of propeller blades all harmonics 

of spanwise loading p „ are antisymmetric about the propeller axis in spanwise 

mN 

coordinate. All, therefore, contribute net moments which drive the wing as 
a whole at the engine station. Were it not for the fact that the frequency 
regime requires that plating response also be considered in detail, these 
loads could be applied in cruder models of wing motion. 

When the number of blades is odd the fundamental distribution is even 
about the propeller axis, contributing a net lift. Higher harmonics are 
alternately odd and even, contributing either torques or forces. 
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D. Predicted Benefit of a Large Number of Blades in Reducing Unsteady 
Wing Airloads Induced by the Wake of an Upstream Propeller 

The argument to demonstrate that for a given value of propeller 
performance thrust and aircraft lift unsteady wing airloads are reduced 
by a large number of blades will here parallel that in Section II. E. 

Eqs. III. 24 and 16 show that at high reduced frequency k (based on 

wing semichord and flight speed) induced wing pressure distributions are 

- 1/2 

proportional ^to purN/R t .k , where purN/R t is in turn proportional to 
T/R with T denoting steady propeller thrust. Since for steady flight the 
amplitude coefficient of wing pressures may be interpreted as proportional 
to the product of the frees tream dynamic pressure and a spanwise- variable 
lift coefficient, we have that for point x,y on the wing surface 


^mN 

steady 

P 


T/R 2 

pU 2 c L <y) 


1 

/k 


(III. 30) 


Moreover, since k=mNftb/U=mN (b/R t ) ( QR t /U) ~mN (b/R t ) ■'•mN , it therefore follows 

that for a given value of design thrust T unsteady wing pressures fall 

- 1/2 

relative to steady pressures as N , as N is increased. 
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APPENDIX A 


JUSTIFICATION FOR EXCLUSION OF MULTIPLE- SCATTERING EFFECTS 


A, 1 INTRODUCTION 

Two of our simplifying assumptions have been that the propeller 
blades effectively lay flat on their plane of rotation without pitch or 
twist, and that unsteady velocities- mutually induced* by unsteady components 
of blade-bound and shed-wake vorticity are negligible in comparison to 
incident modal velocity amplitudes. 

The first assumption causes mutual influences through direct pressure 

fields to be identically zero because flat-plate loadings for each blade 

involve then only thrust dipoles normal to the propeller disk. We shall 

examine its validity by estimating the pressure field with which two 

neighboring blades insonify each other, given that in fact a finite angle 

of twist exists between them. The procedure will be to use the aerodynamic 

theory available for a single airfoil at zero pitch or twist, and then to 

tilt the resulting source distribution through the actual pitch angle to find 

non- zero values of near-field pressure at off-plane points occupied by two 

adjacent blades. The results will be normalized by the input disturbance 

th 

pressure pUw M corresponding to the mN harmonic of the incident velocity 
mN 

field, and their fractional value will represent the pressure field’s measure 
of unsteady aerodynamic interference, or multiple scattering, between blades.* 

Similarly, estimates of interference through induced flows will rely 
on the zero-pitch loading solution appropriately tilted according to pitch. 

The ratio of such results to the amplitude of the incident velocity’s first 
harmonic (say) will indicate a-posteriori to what extent each blade’s 
aerodynamic response can be expected to be independent of others * . 

Fig. A. 1 shows a ’’radial” view of two of the N blades at an arbitrary 
circumferential cross-section through the propeller. As stated above, we 
shall use blade loadings corresponding to passage at zero twist and pitch. 

We shall further simplify the calculation by considering only the leading- 
edge contribution of such loading and its associated pressure field, i.e., 
by taking the airfoil to have its trailing edge at infinity. 

Mutually induced pressures are not used as input conditions in rigid-airfoil 
aerodynamics; however, certain situations may require simple estimates of (cont.) 
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A. 2 BLADE-BLADE AERODYNAMIC INTERFERENCE THROUGH DIRECT PRESSURE FIELDS 

For the semi-infinite chord airfoil Ref. 14 gives the pressure field 
an arbitrary field point (near or far) : 

cos [1/2 tan 1 (6tan0)] 

p(r,0,t) = -puw ■ ■ ■ ' ~ 

/ o 2 2 2 1/2 

fi/irr (p+k^/8 ) (1-M cos A sin 9) 

k cos^A • x 

.exp[i{ oit-k s + — ; - ■ — -Tr/4-ur} ] (A.l) 

S 6 

where in the present application, we have that 

k = k b cosA 
x mN 

a) = S3mN 

k s = mN sin ^ y/R (A* 2) 

s J 

-1 /~2 2 
0 = tan z/x , r =* x +z 

etc, with k . b, and A all functions of radial position R. In Eq. A.l 
mN 2 

U=ftR,M=U/c , 8 = / 1-M cos^A and w T are also functions of R? also, the 9 

o mN 

range must be defined as (0,2rr) , We drop the b-subscript on the blade- 
based coordinate system defined in Fig. II. 3. 

For m=l (the blade-passage fundamental) we have applied Eq. A.l to 

calculate airfoil A's pressure field at airfoil B's location, and vice- 

versa. The results normalized by pw^U are displayed in Fig. A.2 

for three radial stations. The maximum interference pressure is 

shown to be about 16% pUw„ for R=R . for blade A on blade B. This result, 

N tip 

based on an airfoil model which does not satisfy the Kutta condition and 
which in fact yields maximum pressures along the wake direction, should be 
a conservative estimate for the lower values of blade pitch near the blade 
tips. Unfortunately, a correcting trailing-edge solution to prove the above 
assertion by modifying Eq. A.l to satisfy the Kutta condition is not available 
in closed form for arbitrary field points and, at any rate, the value of 16% 


* 

total chordwise gradients as input, for example, to unsteady boundary- layer 
flow calculations. Another possible application is the total loading influence 
function as input in an aeroelastic calculation. 
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is sufficiently low to justify neglecting the forward-scatter from airfoil 

A on airfoil B. The broken lines in Fig, A. 2 summarize the computed levels 

of normalized backscatter (airfoil B on A) . As expected, they are for 

the most part significantly below the solid lines, not just because the 

spurious trailing-edge point is further away, but because compressibility 

effectively lengthens distances to upstream points (as well as enhances 

propagation downstream by washing it along) . The effect is embodied in 

the argument of the cosine factor in Eq. A.l. Let us take, for example, two 

angular positions symmetric about the 8 =tt/ 2 divide between upstream and 

o o 

downstream regions at zero pitch, say 0^=45 for downstream and 0 2 =135 for 

upstream. With g=l, i,e., without compressibility, 0^/2=22.5 and 

0 2 /2=67,5°. With 6=1/2 (M=. 866 for A=0 , however , tan 1 ( 6tan0) becomes Respectively 

for 0^,0 2 # tan ^ 1/2=26.6° so that 1/2 tan 1 (6tan0) =13. 3°, and tan 1 (-1/2) =153. 4° 

so that 1/2 tan 1 ( 6tan0) =76. 7°. Compressibility has therefore increased the 

difference between upstream and downstream signal amplitudes from 

M=0: cos (22.5°) -cos (67.5°)=.924-.383=.54, to M=.866: cos (13. 3°) -cos (76. 7°) 

=.973-. 230=. 74. 

In principle Eqs. A.l and 11.21, etc., hold up to supersonic Mach number 
sec A, when sweep is present. 


It is important to remember that the model of the airfoil with semi- 
infinite chord becomes irrelevant at M=0 for the special case of zero blade 
sweep, and that it is useful in the so-called hyperbolic regime only when 
the frequency and Mach number are both hiqh enough to yield certain minimum 
values of u and k /6^- (cf. Ref. 5's skcwed-gust analysis) . 
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A. 3 BLADE-BLADE AERODYNAMIC INTERFERENCE THROUGH INDUCED FLOWS 


The potential $ associated with the pressure field of Eq , A.l is 
obtained from Bernoulli's equation 


D 4> 

Dt p 


(A, 3) 


where D/Dt is the linearized substantial derivative 3/3t+ (U/b) 3/3x. The 
solution of Eq, A. 3 is conveniently expressed in terms of contour integral 
C (Fig. A. 2a) 


$(x,s,Z,t) 


ibw 


mN 


2tt& 



, exp fi{ u)t-k s s+k^M 2 cos 2 Ax/8^} ] 


. f kl expf -iAx- | Z | A 2 -p 2 } (A>4) 

C (X-k x /B 2 ) 

where in the present context u>,k ,k , etc are as given by Eqs. A. 2, Z is 

x s 

the nondimensional normal to airfoil coordinate compressed by the B factor, 

2 

i.e,, Z=Bz; £ is shorthand for p+k^/8 . 

Differentiating with respect to zb we obtain the ratio of induced 
vertical velocity w 1 {x / s,Z) to incident velocity harmonic at field 

point x,Z. 


w (x,s,Z,t)= - — “ expfi{ojt-k s+k M cos Ax/8 } 1 

2 7T Jc s x 

W K, 

mN 


y* 3X/X + |j~ exp f -iXx- 1 z] A 2 -m 2 } 




If we let Z=0 and take x>0, Eq, A. 5 yields 


w 1 ( x>0 /S yZ-0 / t) 


w 


mN 


: -exp [i{ ojt-k s -k x)l 
s x 


(A, 5) 


(A. 6) 
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as required by the boundary condition of flow tangency on the semi-infinite 
chord extending over CKx<°°. For x<0, Z=0, the branch point contribution 
(Fig. A. 2b) from X=-y gives the vertical velocities induced upstream on 
the plane of the blade: 


-k s) ] 
x s 


wN* o I _ s,z.q x t I ,_ el!p , i(lot _ k 

W mN 

| 1- (l+i) E* [C (-x) ] - exp.[- i(ir/4 + C( -I x 

l /iTt (-X) 


(A. 7) 


For large values of -x, i.e., far upstream, and specializing to the case of 

2 2 

zero sweep so that y= k^M/3 =kM/3 , the above becomes proportional to 


,ikM r . ,, 
exp it— 

1-M 


[5(-x) ] 


3/2 


c t 



indicating a quickly decaying signal of very short wavelength as M+l 

travelling at speed -(l-M)c relative to the airfoil (-c relative to the 

o o 

freestream) . In Eq. A. 7 the vertical induced velocity becomes infinite 
as x-O because at x=0 + the blade-bound vorticity has an infinite value. 

Eq. A. 7 could by itself be used to estimate upstream flows induced 

by blade B on A if one neglects the pitch angle so that the two blades are 

coplanar; small values of I w 1 (x<0,0 # t) /w I over the chord of A would 

mN 

then justify basing flow tangency on incident velocity alone as far as 
upstream influences are concerned. Here, however, we are also interested 
in potentially stronger shed-wake effects, and for this reason will keep 
the two airfoils at nonzero pitch. 

In calculating the shed wake a choice must be made for the model of 
loading on an airfoil of finite chord. Here we choose that given by 
Eq. A. 1 0=0 and corresponding to the leading-edge solution, but applied 
only over 0<x<2. For x>2 the wake condition of p=0=(Ap) will be enforced. 
At sufficiently high frequency or subsonic Mach number, the trailing-edge 
correction to such a loading has the relatively unimportant local effect 
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of "rounding off" the resulting pressure discontinuity at x =2 by providing 
it with the proper / 2-x decay as x-*2 . For simplicity of notation and 
analysis we neglect sweep effects and take k^=k (A=0) . 

The calculation can proceed along the lines of the analysis for the 
upstream influence by considering the singular structure of the velocity's 
/.-transform, or can instead begin by posing Possio's integral for a point 
in the wake. For the sake of variety we take up the latter approach: 

2 

w s (x,z=0,t) = - -77 f ^ ApiC.t)K(|x-5|> (A. 8 ) 

P V J 

0 


with K ( I x- 5 1 ) given by 


K(x-C) = 77 . e 1MM(x C) (iM H (2) [p (x-£) ]-H (2) tu(x-C) 1 ) 

4p i 1 O 


k(x-S)/B 


2 -ik(x-C) 
+ 16 e 


/ 


du (M | u| ) 


(A. 9) 


after normalization by the semichord b and noting that |x-£|=x-£ for 

x>2 (the wake). The superscript s on w S (x) stands for shed wake. In Eq. A. 8 , 

Ap is given by p + -p =2p or twice the result in Eq. A. I with 8=0, 0<r=x<2: 


w (x,z= 0 ,t) 
W mN 


= ( 1— i ) 


h 


2 k 

(1+M)TT 


2 

iut / d£ 

J /r 


-ip ( 1-M) £, 
5 K(x-C) 


(A. 10) 


Incidentally, the upstream result in Eq. A. 7 for a semi-infinite chord 
airfoil can also be obtained from Possio's integral applied for x<0. Then 
we would have that 


x<o : 


w (x,z= 0 ,t) _ 


= ( 1 -i) 


2 k 


loit 


mN 


( 1+M) IT 


oo 

ffi 


-ip (1-M) £ 

e C | ) 


K( x- 


( A . 11) 
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where K(|x~c|) is now given by 


K(5-x) = 


46 


iyM(x-C) , . M u (2) . . (2) ... 

e (- 1 M H fy(£-x)]-H fy(C-x)]) 

1 o 


+ i6 2 e -ik(x-0 


k(x-C)/6 


/ 


du e^ 11 H ^ (-Mu)' 
o 


(A. 12) 


In Eq. A. 11 the integral extends to infinity because the leading-edge 
loading has not been cut off past x=2 as required by Ap(x>2,z=0, t)=0-r-the 
wake condition implicitly applied in Eq. A. 10. 


We begin by looking at points far downstream from the trailing edge 
in terms of wavelengths and later examine in what sense the far wake so 
computed can be expected to hold closer to the trailing edge, given the 
contemplated high reduced frequencies and moderate-to-high subsonic Mach 
numbers. We approximate the Possio kernel by its large kx behavior. Thus, 


k(x-5)/B' 




l 


du H ^ (M } u| ) e 


1U 


= ie 2 e - ik(x -^ 


f-f 


du H (2) (Mu)e 1U 
o 


k (x-£) /6 


- - 26e 


-ik(x-S) 


j>-/f 


U-<i-i)E 'ITS ( *- EI »| 


■26e 


-ik (x-C) + ^-iu(l-M) (x-g) I (1+i) (1+M) 4 o ( ^ _3 /2 

/ Tty (x-£) 


(A. 13) 


Similarly, for y(x-£)>>l we have that 
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(2) l 2) 

iyM(x-S) {iM H l [y(x-C)]-H [y(x-C)]} 

e l o 

-iu(l-M) U-O j (1+i) (1+M) . -3/2., 

~e j- — + 0 (y ) 

' / it y ( x- C ) 

which cancels the second term on the right side of Eq. A. 13 so that 


K(x-£)=-l/2 e -lk ( x -^ + o fy(x-S)}" 3/ ? exp { - ^ (x-£)} (A. 15) 

p ( x- £ ) 


Substituting Eq. A. 15 into A. 10 we finally have that 


kx 


lim w (x,z=0,t) 
w 


mN 


-ik (x-Ut/b) 

-d-i )E (— )e 


(A. 16) 


The exponential factor in Eq. A. 16 carries the obvious physical 

interpretation that the wake is convected and thus steady in the still fluid 

(for which the nondimensional x-coordinate becomes x =x-Ut/b) . As 2k/(l+M)-**> 

s 

(recalling that the 2 stands for nondimensional chord) we have that 


w S (x,z=0,t) 
kx-*° 


e-ik(x-Ut/b) {1+0(k -l/2 )} 


(A. 17) 


This initially surprising result states that at high frequencies the 
"inviscid” far wake shed by an airfoil continues to cancel the vertical 
velocity of the input gust. In a more general high-speed context, it 
suggests the possibility of effectively nulling out for a chosen streamline 
the normal velocity component of an undesirable convected disturbance through 
the introduction of an airfoil. 


There is yet another path to Eq. A. 16 which although perhaps analytically 
less self-contained is at the same time physically more revealing* We write 
the integral equation for the strength Y(* s ) of the vortex sheet in the far 
wake in the frame of reference of the still fluid: 


w S (x ) = - 
s 



dfy ill 


(A. 18) 
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where the integral equation is of the Cauchy-Principal Value type and is 
so designated. The kernel [-2 tt(x s -£)] ^ solves the Laplace equation and is 
the usual special case of Biot-Savart 1 s Law applied to the self-induced 
normal velocity field; it may be obtained from the Possio kernel for zero 
frequency and Mach number. 


Fourier transforming Eq. A. 18, recognizing that 


00 

-i-.i f 
V s 2 J 


hi JL. ~ iX "0 
dX lYf e 


(A. 19) 


and that w(A)=FT{w (>0} is proportional to <5(A-k), we therefore have that 


y ( A) - 4iri 


W 


6 (X-k) 


(A. 20) 


-C .U. ^ 

1.1U1U WliXUi 

y(x ) = 2iw S (x ) => iw s =u + (A. 21) 

s s 

since y=u + -u =2u + , where u*' denote tangential velocities just above and 
and below the shed wake. 

Returning now to the airfoil's frame of reference we recall that 
application of p=Ap=0 along the wake requires that 


DJ> 

Dt 


(x >2 ,o, t) =0 


(A. 22) 


whose solution is 

4>(x,o,t) = <Mx=2,o,t - (A. 23) 

since 0,o, t) =0. With <Mx,o + , t)= $(x,o )exp(iiut) Eq. A. 22, however, 

states that 

iu)$ (x,o) e laJt =- jy(X|t) (A. 24) 

which via Eqs. A. 23 and A. 21 yields 
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y(x,t) = — - — •4>(2,o) expf ioi [t- (x-2) b/U]} 


(A. 25) 


w S (x,t)=- k exp{iw [t- (x-2)b/U]} 


(A. 26) 


Now, in Ref. 5 Eq. 15 Amiet has calculated <t(x,o) to be 


♦<-<» - t b “«N 


(A. 27) 


and so we finally obtain that 


s ov 

w ( x , t ) =“W (l-i)E(-r— — ) expf id) [t-xb/u]} 
mN 1+M 


(A. 28) 


thus reproducing Eq. A. 16, which we now see is the result of two physically 
cancelling effects. On the one hand we have the k 1 loss of lift at high 
frequencies displayed by Eq. A. 27 ($(2,o) is 1/2 the total airfoil bound 
circulation) ; on the other we have the linear dependence on frequency of 
the rate of vortex shedding as shown by Eqs. A. 24-26. 

For very high frequencies far-wake vertical motions have been found 
in Eq. A. 17 to be identical to those self induced by the semi-infinite 
loading of Eq. A. 1. The obvious conclusion is that given a willingness to 
neglect the effects of the near wake, which is reduced in both strength 
and extent with increasing frequency, we may take the solution of Eq. A. 1 
to model the vertical velocity distribution over z=0 for an actual 

airfoil of finite chord. Again, the essential sin committed is omission 
of an additional effective origin introduced by the near wake in the 
vicinity of x=2 + . We justify it by reasoning that next to the infinite 
discontinuity in normal velocity across x=0— shown by Eqs. A. 6, 7 the 
neglected source should be comparatively unimportant. 

Reintroducing sweep, in a form suitable for computation (removing the l//^x 
singularity) , mutually induced normal velocities for airfoils A and B are 
given by 


- — ^ x ■ ' 2 i ^ = exp f i (Tr/2+k M 2 cos 2 Ax/(i"-k s+tot) } pB | z | 
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exp(-i (tt/4-K5 ) 


)} 


(A. 29) 


When calculating the influence of A on B, x=0 is at A's leading edge; for 
that of B on A,x=0 is at B's leading edge. In both calculations |z| is the 
hlade-blade vertical separation given by (2irR/bN) sin 6 t . 

Fig. A. 3 shows predicted induced velocities based on Eq. A. 29. As 
expected they are more pronounced for A on B than vice-versa, with an 
overall high in amplitude ratio of roughly 14%. We conclude that for the 
present situations interference through induced flows is negligible by 
about the same margin as direct pressure fields. 
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(z-z w )/c w X 10 3 


Fig. II. 1 Mean wake measured at 321-chord range downstream of the trailing 

edge of a typical full-scale aircraft. (courtesy of Douglas Aircraft 
Co. ) 
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Fig. II. 2 Idealization for wing interference model (a,b);and coordinate 

system "or propeller blade aerodynamic analysis (cf. Fig. TI. 3 

for expanded view of blade point (T)) . 
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2IIR 

mN 



Fig. II. 3 Definition of geometric parameters at typical blade radial station, 
point(T) in Fig. II. 2c; blade is assumed stationary and gust travels 
in positive 9 direction because blade rotates in negative 9 direction 
(adapted from Ref. 4). 
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A(R) S A( v/R K2 + y 2 ) = tan -1 [ay (R*) /dR*]- tarf 1 [y(R*)/R*] 

4 Schematic of modeled mean-swoop lino for a bl ado, showi ncj 
definitions of sweep and other parameters; global version 
of Fig . IX . 3. 



y/y T = ( R / Rj) 

R*/R* = Rh /r t +o-(1-R h /Rt) 
r h /R? = R h /R t_ 
n/i-CVt/Rt ) 2 



o- y T 


Fig. II. 5 Mathematical model of family of blade-sweep lines, one for 
each value of power n, using y^/R t =l/2, R u /R t =l/4. 
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(b) 

Fig. II. 6 (a) Frontal view of full-scale propfan prototype; 

(b) Expanded view of one of its blades indicating approximate 
circumferential dimensions and circumferential components 
of blade sweep angle. 







Fig. II. 7 Top and side views of Gulfstream IT wing and nacelle showing 
dimensions relevant to calculations. 
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11.14 Predicted spectrum of thrust propeller force (re 1 lb), for propeller 100% of chord 
downstream of pylon trailing edge. 


Fig. 11.15 Predicted inflow fundamental and sectional load vs radial position 
for a wing-mounted pusher propeller; propeller position is 32%-c^ 
downstream of trailing edge; (a) propeller axis on wing plane, 

(b) propeller axis at elevated nacelle position (1.78 ft). 
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at elevated nacelle position (z q -1.78 ft); propeller 





for single vortex 

. H v-vo 

w 2 tt (7-7 0 ) 2 + h 2 
y 0 = R cosilt, h 0 =Rsinilt 



3 | Or 

I 


Fig. III.l Wake model for a four-bladed propeller, indicating horizontal 
plane of wing and vortex-induced downwash. 
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N even 



N odd 



Fig. III. 2 Spanwise behavior of normalized downwash fundamental 4 hR w /TN 
for: (a) even-b laded propeller; (b) odd number of blades. 



Fig. III. 3 Assumed triangular distribution for steady blade lift for 
calculation of vortex strength. 
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Ib/ft 



Fig. III. 4 Predicted pressure distribution 






Ib/ft 



Fig. ill. 5 Predicted pressure distribution for wing ahord 




= 4R 



Fig. III. 7 Predicted pressure distribution for wing chord at y=4R. 



x = l (half chord) 



Fig. III. 9 Predicted spanwise pressure distribution for chord point x=L. 
(1/2 chord) . -74- 





Fraction of blade chord 

Fig. A. 1 Calculated unsteady interference through direct pressure fields, dashed 
lines show backscatter of blade B on A, solid lines show forward scatter 
(blade A on B) . -76- 
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Velocities induced by 
air foi I A on airfoil B 

Velocities induced by 
airfoil B on airfoil A 


.75 1 

Fraction of blade chord 

Fig. A. 3 Calculated unsteady interference through induced flows; dashed lines 
indicate influence of airfoil R on A; solid lines that of A on B. 
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